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Abstract 

Classical solutions of the spherically symmetric Nordstrom- Vlasov sys- 
tem are shown to exist globally in time. The main motivation for investi- 
gating the mathematical properties of the Nordstrom- Vlasov system is its 
relation to the Einstein- Vlasov system. The former is not a physically cor- 
rect model, but it is expected to capture some of the typical features of the 
latter, which constitutes a physically satisfactory, relativistic model but 
is mathematically much more complex. We show that classical solutions 
of the spherically symmetric Nordstrom- Vlasov system exist globally in 
time for compactly supported initial data under the additional condition 
that there is a lower bound on the modulus of the angular momentum of 
the initial particle system. We emphasize that this is not a smallness con- 
dition and that our result holds for arbitrary large initial data satisfying 
this hypothesis. 

1 Introduction 



In astrophysics, systems such as galaxies or globular clusters are often modeled 
as a large ensemble of particles (stars) which interact only by the gravitational 
field which they create collectively. In such systems collisions among the par- 
ticles are sufficiently rare to be neglected. Let f = f(t,x,p) > denote the den- 
sity of the particles in phase-space, where ieR denotes time, igK 3 position, 
and p£R 3 momentum. This density function satisfies the Vlasov equation — a 
continuity equation on phase space — coupled to the field equation(s) for the 
gravitational field. In the non-relativistic case, i.e., for Newtonian gravity, the 
field is governed by Poisson's equation and the resulting system is called the 
Vlasov-Poisson system. The initial value problem for this system is by now 
well understood, and satisfactory global existence results have been obtained, 
cf. |51 11UI mi I17| . In general relativity, the Poisson equation is substituted by 
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Einstein's equations, which coupled to the Vlasov equation yield the Einstein- 
Vlasov system. In contrast to the Vlasov-Poisson system this system is highly 
non-trivial even when matter is left out, since the Einstein equations by them- 
selves constitute a nonlinear system of PDE's. Little is known in general about 
the global structure of vacuum solutions, and even less is known for matter 
spacetimes. However, for small initial data there is a global existence result 
for general, asymptotically flat vacuum spacetimes [51. By imposing symmetry 
conditions on spacetime global existence results are known also for initial data 
unrestricted in size. Some of these results also hold for matter spacetimes when 
matter is described by the Vlasov equation, i.e., for the Einstein- Vlasov system. 
It should be pointed out that similar global results for other phenomenological 
matter models are not known, indicating that in general relativity a kinetic de- 
scription of matter is mathematically convenient. For a review on global results 
for the Einstein equations and for the Einstein- Vlasov system, cf. ^JQ. The 
global existence results that hold for unrestricted data for the Einstein- Vlasov 
system all concern cosmological spacetimes which are spatially compact. In the 
asymptotically flat case, where for example gravitational collapse of an isolated 
body is studied, the global structure of solutions to the Einstein- Vlasov system 
is not known for large initial data, even in the spherically symmetric case. Some 
qualitative information is obtained in On the other hand, for small spheri- 
cally symmetric initial data, spacetime is known to be geodesically complete and 
thus to contain no singularities F° r large data singularities will form — cf. 

— an d the central problem of weak cosmic censorship in general relativity 
can be studied, i.e., to show that formation of singularities will always result in 
black holes. The significance of the latter problem motivates a study of a less 
complex but related model, the Nordstrom- Vlasov system. 



2 The Nordstrom-Vlasov system 



In the present paper we investigate a relativistic model which is obtained by 
coupling the Vlasov equation to the Nordstrom scalar gravitation theory (jj]. In 
this theory, the gravitational effects are mediated by a scalar field cf>, and the 
system reads 



A, 



o 40 



dp 



0tf + p-V:rf- 



(a t 0+p-v^)p+(i+ P 2 )- 1 / 2 v^ 



•V p f = 0. 



(2.1) 
(2.2) 



Here 



denotes the relativistic velocity of a particle with momentum p, p 2 — \p\ 2 , and 
units are chosen such that the mass of each particle, the gravitational constant, 
and the speed of light are all equal to unity. A solution (f , (j>) of this system is in- 
terpreted as follows: The spacetime is a Lorentzian manifold with a conformally 
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flat metric which, in the coordinates (t,x), takes the form 



^ = e 2 *diag(-l, 1,1,1). 

The particle distribution defined on the mass shell in this metric is given by 

f ph (t,x,p) = tft,x,e*p). (2.3) 

More details on the derivation of this model are given in 0. It should be 
emphasized that although the system does not constitute a physically correct 
model it still captures some of the essential features of the Einstein- Vlasov 
system which are not present in the Vlasov-Poisson system. In particular, the 
Nordstrom- Vlasov model does allow for propagation of gravitational waves, even 
in the spherically symmetric case. The hope is that the analysis of this model 
will lead to a better mathematical understanding of a whole class of nonlinear 
partial differential equations and eventually to a better understanding of the 
Einstein- Vlasov system. For previous studies of the Nordstrom- Vlasov system 
we refer to 0JE], where existence of local classical and global weak solutions is 
established, and to [H] where the non-relativistic limit is considered. 

It turns out to be convenient to rewrite the system in terms of the new 
unknowns (f,4>), where / is given by 

f(t,x,p) = e 4 ^(t,x,p). (2.4) 

The system then takes the form 

d*<j>-A x <}>=-p, (2.5) 

»(t,x)= [ f(t,x,p)-=L=, (2.6) 
J ^1+p 2 



Sf- 

where 



(S^p+il+p 2 )- 1 ^^] -V p f = 4fS<p, (2.7) 



S:=d t +p-V 3 



is the free-transport operator. The function fi is the trace of the energy mo- 
mentum tensor. We supply the system with the initial conditions 

f(0,x, P ) = r(x,p), </>(0,x) = 4tf(x), dt4>(Q,x) = <l^(x), x, P eR 3 , (2.8) 

and we assume that the initial data have the regularity 

f m G Cl(R% $ G C 3 (M 3 ), 4 n G C 2 b (R 3 ). (2.9) 

Here the subscript c indicates that the functions under consideration have com- 
pact support, b indicates that they are bounded together with their derivatives 
up to the indicated order. 
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We impose spherical symmetry on the initial data. By uniqueness, spherical 
symmetry propagates from the initial data so that the solution (/,</>) is also 
spherically symmetric, i.e., 

f(t,Ax,Ap) = f(t,x,p), cf>(t,Ax) = <t>(t,x), :r,peK 3 , t>0, AeSO(3). 

By abuse of notation we can then write 

f(t,x,p) = f(t,r,u,a), (f>(t,x)=<f>(t,r) 

where 

r:= u:= \p\, a:=Z(x,p), 
the latter denoting the angle between x and p. In particular, 

x 

V x (j){t,x) = d r <j>{t,r)-. 

r 

Consider the characteristic system 

x=p, p = -{S4>)p-{i+p 2 )- 1/2 v x 4> 

of the Vlasov equation i|2.7[l . Due to spherical symmetry the angular momentum 
is conserved along solutions of the characteristic system. Indeed, defining 

L := r 2 u 2 sin 2 a = \xxp\ 2 , 

we find that along characteristics 

^^.-W) i(a; ( s )^( s ))) = 2e 2 ^(9 t 0+V x 0-p) 

+ 2e 24, {xxp)- pxp-xx((S<p)p+{l+p 2 )- 1/2 d r cl)^ 

= 2e 2 ' t 'LS(j)~2e 2 ' l 'LS(j) = Q. (2.10) 

The initial datum / m is assumed to satisfy the additional condition that there 
is a positive number Lq such that 

L(x,p) = r 2 u 2 sin 2 a>L o >0, (x,p) G supp f in . (2.11) 

In a continuation criterion has been proved for the Nordstrom- Vlasov 

system which says that a classical solution can be extended beyond t = T if the 
support of the momentum can be controlled on the time interval [0,T[. More 
precisely, if the quantity 

P(t) = Bup{|p|:Q<«<t 1 (x,p)esupp/( S )}, (2.12) 

is bounded on the time interval [Q,T[ then the solution can be extended to a 
larger time interval [0,T + <5[ for some S > 0. Hence, if [0,T[ is assumed to be the 
maximal life span of a solution, a bound of P(t) on [0,T[ implies that T = oo. 
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In this work we show that the quantity P(t) can be controlled if the initial data 
satisfy the condition 12.11fl . Our main result can now be formulated. 

Theorem For spherically symmetric initial data (f m ,cf) ln ) which satisfy the 
conditions \2. 9\) and 12.11]) . there exists a unique global classical solution 
(/,0)eC 1 ([O,oo[xR 6 )xC 2 ([O,oo[xR 3 ) of the Nordstrdm-Vlasov system fE3|) . 
$2.6)) . \2. 7[ ) ; which satisfies the initial conditions }2.ti\) . 

For the proof it will be essential to look at the behaviour of another quantity 
along characteristics (besides the angular momentum): Along characteristics, 

A(eVlT7)=^^=L=. (2.13) 
ds\ J \/l+p 2 

It is of interest to note that this relation is independent of the symmetry as- 
sumption; the quantity under consideration is the particle energy which would 
be conserved for a static solution. 



3 Proof of the theorem 

We assume that the length of the maximal existence interval T is finite and 
show in three steps that the quantity P from l|2.12|l is bounded on [0,T[ which 
then contradicts the continuation criterion joj Prop. 2 ]. 
Stepl: An estimate for dt<f> 
We split (f> as follows: 

0=0hom+V'- (3.1) 

Here 0hom is the solution of the homogeneous wave equation with the given 
data, and ip is the solution of the inhomogeneous wave equation with zero data. 
Clearly, 

||0homW||oo + ||a t hom (i)|| oo <C, 0<t<T; (3.2) 
recall that we assume T < oo . By 7 , Lemma 7] , 

f dy 1 /** r r +t~ T 

ip(t,x) = - n(t-\x-y\,y)- - = -— / / (i(T,\)\d\dT. 

J\x-y\<t F-J/l Wo J\r-t+r\ 

Note that < so that 4> < </>hom < C. Distinguishing the cases < r < t and r>t 
we get 
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dt~>P{t,r)<— I n(T,t-r-T)(t-r-T)dr, (3.3) 



i) 



where (t — r)+ denotes the positive part of t — r. Define 

P(t) := sup{e^^ |p| \(x,p) G supp/( S ), < s < t} 
Then by conservation of angular momentum, 

/>0O />7T 

rfi(t,r) = 2tt / f(t, r,u,a)r sinauda . 

Jo Jo Vl + w 2 



du 
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OO /»7T 



< Ce~^ r) / / f(t,r,u,a)dadu 

Jo Jo 

= Ce 3<Kt,r) / e -**(*.»-)/(t,r,u,a)dadu 

■/o ./o 

< Ce 2 ^'' r )p(t)<(7P(t); (3.4) 
recall that < </>hom < C and 

f(t,x,p) = f ia (X(0),P(0))exp[4cf>(t,x)-4(^(X(0))] <C, (3.5) 

where (X(s),P(s)) are the characteristics of the Vlasov equation with 
(X,P)(t) = (x,p). If we substitute (|3.4|l into l|3.3[) we get the estimate 

dt<if>(t,r)<-P(t), r>0, 0<t<T. (3.6) 
r 

Step 2: A bound on P 

With and i|31)|> equation (|2~T3|) implies that 



— (e^^/l+p 2 ) <C + CP(s)—< C + CP(s) 
ds V / to 



e 0( s . r )r'usinQ; 



Here we used the fact that < C . At this stage we recall the assumption (|2.11|l 
on the data, namely 

L(x,p) = r 2 u 2 sin 2 a > L > 0, (x,p) G supp/ in . (3.7) 

Using the fact that e 2 ^L is constant along characteristics we get 

^-(e^^/TTp 2 ~)<C + CP( S ), 

and by Gronwall, P is bounded. 
Step 3: Completion of the proof 
Using the bound on P and (|3.5|l we can bound /i, 

H(t,x)= [ f{t,x,p)^L=dp<Ce A ^ j dp<C, 

J e Ht,*)\ p \<c y/l+p* Je*(*-*3|p|<C 



and hence also 0, 



|0(t,a;)| <C+C I r^< C - 



x-y 



<t \x-y\ 



Hence the bound on P implies a bound on \p\ over the support of /, by the 
continuation criterion [SJ Prop. 2] global existence follows, and the proof of the 
theorem is complete. 
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